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Mn doped semiconductors are extremely interesting systems due to their novel magnetic properties 
suitable for the spintronics applications. It has been shown recently by both theory and experiment 
that Mn doped GaN systems have a very high Curie temperature compared to that of Mn doped 
GaAs systems. To understand the electronic and magnetic properties, we have studied Mn doped 
GaN system in detail by a first principles plane wave method. We show here the effect of varying 
Mn concentration on the electronic and magnetic properties. For dilute Mn concentration, d states 
of Mn form an impurity band completely separated from the valence band states of the host GaN. 
This is in contrast to the Mn doped GaAs system where Mn d states in the gap lie very close to 
the valence band edge and hybridizes strongly with the delocalized valence band states. To study 
the effects of electron correlation, LSDA+U calculations have been performed. Calculated exchange 
interaction in (Mn,Ga)N is short ranged in contrary to that in (Mn,Ga)As where the strength 
of the ferromagnetic coupling between Mn spins is not decreased substantially for large Mn-Mn 
separation. Also, the exchange interactions are anisotropic in different crystallographic directions 
due to the presence or absence of connectivity between Mn atoms through As bonds. 

PACS numbers: 75.50.Pp,75.70.-i,71.70.Gm 



I. INTRODUCTION 

Diluted magnetic semiconductors (DMS) are consid- 
ered to be potential candidates for present and future 
technological applications in semiconductor spintronics 
0. During the last decade, there have been numerous 
experimental and theoretical studies of II- VI, III-V and 
IV- VI DMS. Among the III-V DMS, Mn doped GaAs 
system has been studied rigorously for the last few years. 
This system shows a Curie temperature (Tc) of 110 K 
for a Mn doping concentration of 10 %. Despite several 
attempts, Tc couldn't be raised beyond 175 K. Recently, 
there have been reports of some room temperature DMS. 
They include Mn doped GaP Mn doped chalcopyrite 
CdGePa 1, Mn doped GaN 0, ||] etc. The origin of 
ferromagnetism in these compounds is still under debate 
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Dietl et aL0 predicted theoretically a high Curie tem- 
perature (- 400 K) for Mn doped GaN (5 % Mn). Their 
theory was based on a mean field model of hole medi- 
ated ferromagnetism. As the Curie temperature for a 
Mn doped GaAs system is comparatively lower, the the- 
oretical prediction for a higher Tc drew much attention. 
Sasaki et al. grew wurtzite Mn doped GaN films by 
the molecular beam epitaxy method. Magnetic measure- 
ments showed a very high Curie temperature of about 
940 K. They ruled out the possibility of phase segrega- 
tion of some ferromagnetic compound e.g. MnGa and 
Mn4N which also have high Curie temperatures. Room 
temperature ferromagnetism in Mn doped GaN was also 
observed by Reed et al. |5i]. Deep level optical spec- 
troscopy measurements show that Mn forms a deep 
acceptor level at 1.42 eV above the valence band maxi- 
mum for GaN doped with small concentration of Mn. It 
is to be noted that Mn forms an acceptor level at 0.11 
eV above the valence band maximum in the case of Mn 



doped GaAs. So, the overlap of the Mn d-states with the 
valence band is rather strong in Mn doped GaAs com- 
pared to Mn doped GaN. 

Recently, there have been a few first principles elec- 
tronic structure calculations of Mn doped GaN systems. 
Fong et al. |2| performed electronic structure calculations 
of Fe and Mn doped GaN using the tight-binding lin- 
earized muffin tin orbital (TB-LMTO) method. Sato and 
Katayama-Yoshida Ji| performed KKR-CPA (Korringa- 
Kohn-Rostoker-Coherent Potential Approximation) cal- 
culations to study the relative stabilities of ferromag- 
netic and spin glass phases. They showed that for a low 
concentration of Mn, ferromagnetism is favored whereas 
for the high concentration, the spin glass phase is sta- 
ble. The disordered local moment model was assumed 
to describe the spin glass phase. They explained the 
origin of ferromagnetism in these systems by a competi- 
tion between double exchange and superexchange inter- 
actions. Kulatov et al. [Tll | studied electronic, magnetic 
and optical properties of zinc-blende (Mn,Ga)N for dif- 
ferent concentrations of Mn by the TB-LMTO method in 
a supercell approach. Anomalous exchange interactions 
in III-V DMS were found from calculations by Schilf- 
gaarde and Mryasov |0 . They predicted aggregation of 
ma gnet ic nanoclusters inside the III-V host. Kronik et 
al. [13 considered (Mn,Ga)N in wurtzite structure and 
performed electronic structure calculations using a plane 
wave pseudopotential method. In a recent preprint jl4| . 
the different origins of ferromagnetism in (Mn,Ga)As and 
(Mn,Ga)N systems were discussed. The authors pointed 
out from self-interaction corrected (SIC) pseudopoten- 
tial calculations that (Mn,Ga)N is characterized by lo- 
calized Mn 3c? states with a strong self-interaction. In 
(Mn,Ga)As, d states are weakly correlated and are rather 
delocalized being strongly hybridized with As p states. In 
this communication, we attempt to understand the elec- 
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tronic structure and magnetic interactions in Mn doped 
GaN and GaAs systems. The motivation of this paper 
is twofold : (a) to investigate electronic structure and 
magnetism of Mn doped GaN system in detail and (b) 
to have a comparison with Mn doped GaAs system. The 
paper is organized as follows : In the next section, we de- 
scribe the computational details. In the results section, 
a subsection describes the electronic structure and mag- 
netism of the (Mn,Ga)N with varying Mn concentration 
in the local spin density approximation (LSDA) . Then we 
present results from LSDA-t-U calculations. Finally, we 
show a comparison of exchange interactions in (Mn,Ga)N 
and (Mn,Ga)As systems. 



II. COMPUTATIONAL DETAILS 

GaN can be grown both in the zinc-blende and wurtzite 
structures. But, usually, Mn is doped in a wurtzite GaN 
host 0]. In our calculations, we have considered the 
wurtzite structure. Also, for a comparison, we show cal- 
culations for the zinc-blende structure. Experimental lat- 
tice parameters such as a=3.189 A and c=5.185 A with 
a c/a ratio of 1.626 were taken for the calculations in the 
wurtzite structure. Results from the atomic relaxations 
revealed that the nearest neighbor bond-lengths between 
Mn and N change only by 3 % compared to that of the 
bulk GaN. This is in agreement with the results of Kronik 
et al. |13]- In general, Mn doping in substitutional site 
results in a small relaxation [l^ of the nearest neighbor 
anions around it. 

Calculations have been performed by an ab-initio plane 
wave code (VAST) Vanderbilt 0| type ultrasoft 

pseudopotentials were used for the LSDA calculations. 
LSDA-I-U calculations were done in the projector aug- 
mented wave (PAW) [l^ method as implemented by Ben- 
gone et al. [13 . Ga 3d orbitals were included in the basis 
set of the Ga pseudopotential. A kinetic energy cut-off 
of 350 eV was used for the plane waves included in the 
basis set. Ceperley and Alder exchange-correlation 
functional parameterized by Perdew and Zunger |2lj | was 
considered within LSDA. We have also checked that the 
results obtained within GGA (generalized gradient ap- 
proximation) p3 | are similar. A 8x8x6 k-p oints grid was 
used in the Monkhorst Pack scheme [23 for small su- 
percells. For the largest supercell considered, a 2x2x1 
grid was used. Local properties such as local density of 
states and local magnetic moments were calculated by 
projecting the wave functions onto spherical harmonics 
|24l |. The radii chosen for the projection were 1.31, 1.21 
and 0.74 A for Mn, Ga and N respectively. 

We have modeled the system using different su- 
percell sizes to simulate different Mn concentrations. 
For wurtzite structure, supercells having 4, 8, 16, 32, 
72 and 108 atoms were used to model a composi- 
tion Mn^Gai_^N for a;=0.5, 0.25, 0.125, 0.0625, 0.028 
and 0.018 respectively. For zinc-blende structure, we 
used a 64 atom-cell to simulate a Mn concentration of 
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FIG. 1: Spin resolved density of states of (MniGai_a;)N in 
the wurtzite structure. Data are shown for (a) pure GaN (b) 
6.25 % Mn in GaN (c) 12.5 % Mn (d) 25 % Mn (e) 50 % Mn. 
In (a), total DOS/cell as well as the p-DOS of N have been 
shown whereas in (b) to (e), d-DOS of Mn (in shade), p-DOS 
of N and total DOS have been plotted. Energies are plotted 
with reference to valence band maximum (VBM) of GaN in 
(a) and Fermi energies {Ef) in (b) to (e). 



x=0. 03125. It has been found experimentally that Mn 
occupies the Ga-site Therefore, in our calculational 
unit cell, 1 Ga atom was substituted by a Mn atom. 



III. RESULTS 
A. Calculations within LSDA 

In a III-V semiconductor, a cation vacancy creates 3 
holes in the valence band leaving anion dangling bonds. 
When Mn occupies the cation site, it donates 3 electrons 
to fulfill the bonding. Mn is left with 4 unpaired d elec- 
trons which give rise to 4 /.t^/Mn atom. In a realistic 
situation, there can be compensating donors e.g. As an- 
tisites and interstitial Mn atoms {2^], present in the sys- 
tem to increase or decrease the magnetic moments. But 
as we are dealing with ideal uncompensated systems, we 
always obtain a magnetic moment of 4 /i^/Mn atom for 
substitutional Mn. 

In Fig. l(a-e), we show the density of states (DOS) 
for various concentrations(a;) of Mn in Mn2,Gai_a;N in 
the wurtzite structure. In (a), the DOS of the undoped 
GaN is presented. The large band gap is evident from 
the figure. The calculated direct band gap at the F point 
is 1.9 eV which is underestimated compared to the ex- 
perimental band gap of 3.4 eV. This well known un- 
derestimation is inherent in the formulation of density 
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TABLE I: Projected charges and magnetic moments (in /is). 
Mn-d'' indicates charge in the the sphere around Mn for d 
electrons. Mrimom and Nmom indicate projected magnetic 
moments inside Mn and N spheres respectively. 



X in Mna;Gai_a;N 


Mn-d'«' 




4 X^nri Nmom 


0.018 


5.04 


3.40 


0.016 


0.028 


5.04 


3.41 


0.016 


0.0625 


5.05 


3.43 


0.01 


0.125 


5.05 


3.42 


0.015 


0.25 


5.04 


3.47 


0.02 


0.50 


5.09 


2.61 


-0.02 



functional theory and is well documented in existing lit- 
erature. From Fig. 1(b) to 1(e), we show the DOSs with 
increasing x. The Mn impurity d peak in the band gap 
is away from the valence band for small x. As the con- 
centration increases, this peak is broadened due to the 
overlap of Mn d wavefunctions and the gap between the 
impurity peak and valence band edge vanishes. The to- 
tal DOS/cell gradually regains the shape of the DOS of 
undoped system as x is reduced. However, up to x=0.25, 
the Fermi level cuts only the spin up DOSs. As there 
is no state at Fermi level for the spin down channel, we 
obtain a half metallic solution giving rise to an integer 
magnetic moment of 4 /i^/Mn atom. For x—0.5 shown 
in Fig. 1(e), the system behaves like a ferromagnetic 
metal with a reduced magnetic moment. One notice- 
able difference between Mn doped GaAs and Mn doped 
GaN systems is the position of Mn impurity band in the 
gap. In (Mn,Ga)As, Mn impurity band almost merges 
with the top of the valence band (approximately 0.1 eV 
above the valence band edge) whereas in (Mn,Ga)N, the 
Mn impurity band is separated from the valence band by 
0.56 eV. The width of this impurity band is 0.94 eV. 
As the spin polarized Mn impurity band is distinctly 
in the gap, the valence band is less spin polarized in 
(Mn,Ga)N compared to (Mn,Ga)As. Spin polarization 

i^VBM/CBM - ^vbm/cbm) valcncc and conduc- 

tion bands at the F point is 0.05 eV and -0.42 eV respec- 
tively for (Mno. 0625 Gao. 9375 )N. e is the eigenvalue and 
VBM and CBM represent valence band maximum and 
conduction band minimum respectively. 
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FIG. 2: Spin resolved density of states of (Mn,Ga)N for (top) 
wurtzite and (b) zinc-blende structures. Here, the Mn con- 
centration is 6.25 %. 

neighbor N atoms are also tabulated. In most of the 
cases, the moments are parallel to Mn moments which is 
not the case for a Mn doped GaAs system where nearest 
neighbor As moments are antiferromagnetically coupled 
to the Mn magnetic moment. 

As GaN can be grown in both zinc-blende and wurtzite 
structures, we have also done calculations for a 6.5 % Mn 
doped GaN system in zinc-blende structure. In Fig. 2, 
DOSs for both structures are shown. DOS for a zinc- 
blende structure is similar to that calculated by Kulatov 
et al. Also the magnetic moment on Mn atom (3.4 

/is) agrees very well. There is no striking difference in 
the broad features of the DOSs for the two structures. 
In the zinc-blende structure, the peak at the Fermi level 
is sharper and the nature of these states are different (t2 
compared to e states for the wurtzite structure). See 
fig. 5 and related discussions. For both the structures, 
the Mn impurity peak in the energy gap of the host is 
separated from the valence band of the host. 



In Table I, we present the local charges and magnetic 
moments of Mn and N atoms for different concentrations 
of Mn. Charges and magnetic moments of Mn remain 
almost the same with the concentration variation of Mn. 
The insensitivity of magnetic moment with concentra- 
tion is a signature of localized d-states of Mn. The total 
moment/cell is always 4.0 /is which is the signature of 
a half-metallic solution. The only exception is the case 
of a;=0.5, where the total moment/cell is 2.77 hb- The 
exchange splitting of Mn d-states is less in this case al- 
lowing both spin-up and spin-down d-states to cross the 
Fermi level. The averaged induced moments on nearest 



B. LSDA+U calculations 

It is a matter of debate whether the itinerant band 
model or the localized atomic model is appropriate for 
the description of Mn doped semiconductors. Density 
functional calculations Esj based on LSDA or 

GGA provide the basis of itinerant picture whereas oth- 
ers models are based on a localized atomic picture . A 
recent photoemission experiment |30| on Mn doped GaAs 
system revealed the main Mn d-peak to be situated 3.4 
eV below the Fermi level. A previous photoemission ex- 
periment [U reported the peak to be 4.4 eV below the 
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TABLE II: Projected charges and magnetic moments (in /is) 
calculated within LSDA+U. The notations are same as in 
Table I. 





U=4.0 eV 


U=5.0 eV 


U=6.0 eV 


U=7.0 eV 


Mn-d*^ 


5.01 


5.00 


4.99 


4.98 




3.72 


3.83 


3.93 


4.02 




-0.02 


-0.04 


-0.05 


-0.06 



Fermi level. Also, the importance of electron correla- 
tion effects in these systems was highlighted. Park et. 
al. [sj argued from LSDA-I-U calculations that correla- 
tion corrections are important to have a better agreement 
with photoemission spectra. On the other hand, all den- 
sity functional calculations based on LSDA show a peak 
around 2.6-2.9 eV below the Fermi level |ll|23|. So it 
can be argued that the completely localized picture or 
the completely itinerant picture cannot solely describe 
these systems satisfactorily. 



We have done LSDA-I-U calculations to investigate the 
effect of electron correlations. Firstly, these calculations 
do not exist in literature and secondly, it is interesting 
to compare this simple technique with more rigorous SIC 
calculations. In Fig. 3(a-d), we show the DOSs obtained 
from LSDA and LSDA+U calculations for 6.25 % Mn. 
As the value of U for Mn d states is not obtained self- 
consistently from first principles calculations, we varied 
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FIG. 3: (a) and (b), Spin resolved Mn d and N p DOSs 
from the LSDA-I-U calculations for different values of U for 
(Mn,Ga)N system in the wurtzite structure, (c) Mn d-DOS 
from both LSDA and LSDA-I-U calculations for (Mn,Ga)N in 
the wurtzite structure, (d) same as (c), but for (Mn,Ga)As in 
the zinc-blende structure. For all plots, Mn concentration is 
6.25 %. 



U from 4 eV to 7 eV treating it as a parameter. In all 
cases, the exchange parameter J was considered to be 1.0 
eV. Increasing U results in a slow shift of the spin-up im- 
purity band towards the valence band. The small peak 
around 1.5 eV below the Fermi level [Ep) is diminished 
gradually in magnitude whereas the peak around -4.5 eV 
below Ep increases in magnitude. Even for U=10 eV (not 
shown here), the localized peak remains pinned close to 
Ep and is not merged with the delocalized valence band 
states. The spin down DOS shifts almost rigidly away 
from the Fermi level towards higher energy with increas- 
ing U. Local charges and magnetic moments of Mn and 
nearest neighbor N atoms are listed in Table II. Mag- 
netic moment of Mn increases with U due to increase 
in localization of d states. The induced moment of N 
atoms also increase. In Fig. 3(c), a comparison between 
LSDA and LSDA+U calculations is shown. The redis- 
tribution of weights of Mn d peaks with the inclusion of 
U is visible. In LSDA+U result, the sharp peak is very 
close to the Fermi level in the impurity band whereas 
it is in the valence band in the LSDA calculation. This 
peak is of e-character and doesn't take part in the bond- 
ing with the neighboring N atoms. These results are in 
agreement with the more sophisticated SIC calculations 
reported by Fihppetti et al. They also found a 

flat band with character at Ep. The nature of the 
states close to Ep from our calculation is shown later in 
fig. 5. To our knowledge, the valence band photoemis- 
sion spectra of (Mn,Ga)N is not available in the litera- 
ture. So, the extent of validity of LSDA approach can- 
not be tested. Our LSDA+U findings can be compared 
with future angle-resolved photoemission experiments on 
(Mn,Ga)N to verify the existence of the localized peak 
close to Ep. A comparative study with Mn doped GaAs 
is shown in Fig. 3(d). The main broad Mn d peak around 
2.8 eV below the Fermi level in an LSDA calculation is 
shifted 4 eV below the Fermi level with a smaller band 
width. The DOS at the Fermi level is also decreased com- 
pared to that of an LSDA calculation. In both LSDA and 
LSDA+U calculations, the hybridization between Mn d 
and As p states are seen. So the holes in the valence band 
have hybridized p — d character. 

In Fig. 4, we show the integrated magnetization den- 
sity around various atoms in the unit cell as a function 
of the radius of integration. It has been calculated as 



M{R) 



[p^ {r) — p'^(r))dr 



where (r) and (r) are the spin- up and spin-down 
charge densities respectively and M is the magnetic 
moment obtained for a radius R. For (Mn,Ga)N and 
(Mn,Ga)As in the zinc-blende structure, the magnetic 
moment reaches the value 4 /i^ at a distance of 5 A far 
from the Mn center. But, in the wurtzite structure of 
(Mn,Ga)N, this value is reached at a smaller distance. 
This again shows a more localized character of Mn d 
states in wurtzite GaN. The integration around nearest 
neighbor anions reveal that N has a positive contribution 
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FIG. 4: Magnetic moment vs. radius of integration around 
specific atoms. In (a) integrated magnetization density 
around Mn spfieres and in (b) same around nearest neigli- 
bor anions and next nearest neigfibor Ga atoms are shown. 
See text for details. Tlie fiorizontal line in the left figure in- 
dicates a magnetic moment of 4 hb- ZB and W represent 
zinc-blende and wurtzite structures. 
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FIG. 6: Total energy difference AE between ferromagnetic 
and antiferromagnetic configurations vs. Mn-Mn distance in 
the unit cell. Data for zinc-blende(ZB) GaN, ZB GaAs and 
wurtzite (W) GaN are shown. The filled triangle represents 
AE for [001] direction in ZB GaAs. 



in (Mn,Ga)N whereas, As has a negative contribution. 
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FIG. 5: Character of states within an energy interval Ef 
to Ef-0.5 eV for (left) (Mn,Ga)N in the wurtzite structure, 
(middle) (Mn,Ga)As in the zinc-blende structure and (right) 
(Mn,Ga)N in the zinc-blende structure having 6.25 % Mn 
in all cases. In x-axis, the plotted states are in the order 
from left : a,p^,py,pz,dxy,dyz,dzx,d^2,d^2_y2. See text for 
details. For the middle panel, both spin up and spin down 
contributions are shown whereas for the left and right panels, 
there is no state in the spin down channel. In the left column, 
the empty bars represent LSDA+U values. For each panel, 
data are shown for Mn(top), four nearest neighbor (nn) anions 
(middle) and the collective contributions from all other atoms 
(bottom) in the unit cell. 

In Fig. 5, we show the character of the states within an 
energy interval close to the Fermi level. The states close 
to Ep in the electronic structure are important in charac- 
terizing the origin of ferromagnetism. The projection of 
the wavefunction onto spherical harmonics around each 



atom has been calculated as mentioned in ref. |2J| . In the 

figure, X is defined as : xfrn = Enfc l(y/m''ICfc)P. where 
a and N are the spin and atom indices respectively and 
n is the index for the bands within the specified energy 
interval. It is clear from the figure that for (Mn,Ga)N 
in the wurtzite structure (left panel) , the d^^ component 
of Mn e-orbital and the pz component of nearest neigh- 
bor N p orbital are the dominant states close to Ep- On 
the other hand, Mn t2 and nearest neighbor As p states 
hybridize strongly in (Mn,Ga)As. This is also true for 
(Mn,Ga)N in the zinc-blende structure. Symmetry of 
the crystal structure and the splitting of d-states under 
the corresponding crystal fields determine the position of 
t2 and e states. For the rest of the atoms in the unit cell, 
states having s and p character are the dominant. 



IV. INTERATOMIC EXCHANGE 
INTERACTIONS 

To determine the interatomic exchange interactions, 
we followed a simple model. In the unit cell, two Mn 
atoms were placed in various positions. For each Mn- 
Mn separation, ferromagnetic (FM) and antiferromag- 
netic (AFM) alignments of Mn spins were considered. 
The total energy difference A£; {AE=EfJ'^-EfJ^) be- 
tween these two alignments is a measure of interatomic 
exchange interaction. In Fig. 6, we plot AE as a func- 
tion of Mn-Mn separation d for both Mn doped GaAs 
and Mn doped GaN systems. For Mn doped GaN, both 
zinc-blende and wurtzite structures were considered for 
these calculations with two Mn atoms in unit cells hav- 
ing 64 atoms and 72 atoms respectively. For Mn doped 
GaAs, a 64 atom unit cell was chosen. Ferromagnetic 
interaction between Mn spins is favored for all the cases 
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considered here. In (Mn,Ga)N, the first nearest neighbor 
(nn) exchange interaction is the strongest. The value of 
IS.E is increased a little bit compared to the same for 
(Mn,Ga)As. So, for a defect-free calculation, there is 
no indication that (Mn,Ga)N should have a much higher 
Tc than (Mn,Ga)As. It suggests that the formation of 
other phases during the growth is responsible for very 
high Tc observed in certain experiments. It is also seen 
from the figure that /S.E decreases sharply with d for 
Mn doped GaN systems. It shows that ferromagnetic 
exchange interaction in (Mn,Ga)N is short-ranged. This 
indicates that the formation of Mn clusters within a short 
radial distance might lead to high values of Tc- For Mn 
doped GaAs, the exchange interaction is long ranged and 
doesn't decrease rapidly. The proper range of ferromag- 
netic interactions can be studied with a bigger supercell. 
In summary, the results for exchange interactions indi- 
cate that the ferromagnetic interaction between Mn spins 
in (Mn,Ga)As is mediated by delocalized valence band 
holes whereas the origin of ferromagnetism in (MmGa)N 
may result from a double-exchange mechanism [l(j| in- 
volving the hopping of Mn d electrons. The other proba- 
ble mechanism can be the formation of Zhang-Rice mag- 
netic polaron (33j . 

Another interesting observation is the anisotropy of the 
exchange interactions in different crystallographic direc- 
tions. In (Mn,Ga)As, ferromagnetic coupling is stronger 
either in the bonding direction e.g. [Ill] or in a direction 
where the two Mn spins are connected by As bonds e.g. 
in [110] direction. The coupling decreases for the [001] 
direction where there is no As atom in between (shown 
as a filled triangle in fig. 6). The strong ferromagnetic 
interactions in [110] or [111] directions result from strong 
p—d hybridization between Mn d and As p states. As the 
interaction is mediated by the delocalized states, they are 
sufficiently long-ranged. In [001] direction, the exchange 
interaction between the Mn spins can not be mediated by 
the p — d hybridization. So, the value of AE is decreased 
compared to the other cases. In (Mn,Ga)N, the exchange 
interaction is not mediated by the delocalized valence 
band states and the ferromagnetic interactions decrease 
sharply with the Mn-Mn separation. So, the anisotropy 
in the exchange interactions is not significantly observed. 
We can conclude that in general, this anisotropy should 
be present for all Mn doped semiconductors where ferro- 
magnetic long range interaction is mediated by delocal- 
ized valence band states. We also studied the exchange 



interactions in (Mn,Ga)N within the LSDA-I-U scheme 
in a similar way described above. U=4 eV and U=7 eV 
were considered for the calculations. We found that /S.E 
for nearest neighbor Mn-Mn distance increases from 170 
mcV/Mn to 240 meV/Mn while going from U=4 eV to 
U=7 cV. Ferromagnetic interactions between Mn spins 
become stronger as the localization of the Mn d-states is 
increased. 

V. CONCLUSION 

We have studied the electronic structure and mag- 
netism of Mn doped GaN systems for a wide concen- 
tration range of Mn. The deep acceptor level of Mn lies 
distinctly in the gap of GaN, separated from the valence 
band of host GaN. This is in contrast to the case of 
(Mn,Ga)As where the Mn forms shallow acceptor level 
close to the valence band of GaAs. Ferromagnetic inter- 
actions are short-ranged in (Mn,Ga)N systems whereas 
they have delocalized itinerant character for (Mn,Ga)As. 
Also, in (Mn,Ga)As, exchange interactions are signifi- 
cantly anisotropic in different crystallographic directions. 
This has been explained in terms of anisotropic p — d 
hybridization. The results presented here are for ideal 
systems having no contribution from the native defects 
formed during the non-equilibrium growth. To have a 
more realistic picture, one should take into account these 
effects. Also, Mn can occupy interstitial positions in the 
lattice and the spin interactions can become quite com- 
plicated. We are presently calculating the formation en- 
ergies of the defects and subsequently the magnetic in- 
teractions in presence of them. Results will be reported 
in future communications. 



Acknowledgments 

We are grateful to the Swedish Foundation for Strate- 
gic Research (SSF) , the Swedish Natural Science Re- 
search Council (NFR), and the Goran Gustafsson Foun- 
dation for financial support. We thank J. Hafner and G. 
Kresse for giving us the possibility to use VASP. B.S. is 
grateful to Dr. Lars Nordstrom and Prof. OUe Eriksson 
for valuable suggestions. Also, we acknowledge compu- 
tational support from National Supercomputer Center 
(NSC) at Linkoping, Sweden. 



[1] Semiconductor spintronics and quantum computation, 

edited by D. D. Awschalom, D. Loss and N. Samarth, 

Springer, Berlin (2002). 
[2] N. Theodoropoulou, A. F. Hebard, M. E. Overberg, C. 

R. Abernathy, S. J. Pearton, S. N. G. Chu, and R. G. 

Wilson, Phys. Rev. Lett. 89, 107203 (2002). 
[3] G. A. Medvedkin, T. Ishibasi, T. Nishi, K. Hayata, Y. 

Hasegawa and K. Sato, Jpn. J. Appl. Phys. 39, L949 



(2000). 

[4] T. Sasaki, S. Sonada, Y. Yamamoto, K. Suga, S. Shimizu, 
K. Kindo and H. Hori, J. Appl. Phys. 91, 7911 (2002). 

[5] M. L. Reed, N. A. El-Masry, H. H. Stadelmaier, M. K. 
Ritums, M. J. Reed, C. A. Parker, J. C. Roberts and S. 
M. Bedair, Appl. Phys. Lett. 79, 3473 (2001). 

[6] P. Kacman, Semicond. Sci. Technol. 16, R25 (2001). 

[7] T. Dietl, H. Ohno and F. Matsukara, Phys. Rev. B 63, 



7 



195205 (2001). 

[8] R. Y. Korotkov, J. M. Gregie and B. W. Wessels, Appl. 

Phys. Lett. 80, 1731 (2002). 
[9] C. Y. Fong, V. A. Gubanov and C. Boekema, J. Elec- 
tronic Materials, 29, 1067 (2000). 

[10] K. Sato and H. Katayama-Yoshida, Semicond. Sci. Tech- 
no!. 17, 367 (2002). 

[11] E. Kulatov, H. Nakayama, H. Mariette, H. Ohta and Yu. 
A. Uspenskii, Phys. Rev. B 66, 045203 (2002). 

[12] M. V. Schilfgaarde and O. N. Myrasov, Phys. Rev. B 63, 
233205 (2001). 

[13] L. Kronik, M. Jain and R. Chehkowsky, Phys. Rev. B 
66, 041203 (2002). 

[14] A. Fi Uppetti, N. A. Spaldin and S. San- 
vito, |cond-niat/0302178| (2003), avahable at 
http : / / xxx.lanl.gov" 

[15] S. Mirbt, B. Sanyal and P. Mohn, J. Phys. Condens. Mat- 
ter. 14, 3295 (2002). 

[16] G. Kresse and J. Hafner, Phys. Rev. B 47, RC558 (1993); 
G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

[17] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

[18] P.E. Blochl, Phys. Rev. B 50, 17953 (1994). 

[19] O. Bengone, M. Alouani, P. E. Blochl and J. Hugel, Phys. 

Rev. B 62, 16392 (2000). 
[20] D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 556 

(1980). 

[21] J. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 



[22] J. P. Perdew and Y. Yang, Phys. Rev. B 45, 13244 (1992). 
[23] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 
(1976). 

[24] A. Eichler, J. Hafner, J. Furthmiiller and G. Kresse, Sur- 
face Science 346, 300 (1996). 

[25] M. Sato, H. Tanida, K. Kato, T. Sasaki, Y. Yamamoto, 
S. Sonoda, S. Shimizu and H. Hori, Jpn. J. Appl. Phys. 
41, 4513 (2002). 

[26] L. Bergqvist, P. A. Korzhavyi, B. Sanyal, S. Mirbt, L 
A. Abrikosov, L. Nordstrom, E. A. Smirnova, P. Mohn, 
P. Svedlindh and O. Eriksson, Phys. Rev. B 67, 205201 
(2003). 

[27] S. Sanvito, G. Theurich and N. A. Hill, J. Superconduc- 

tivity:Incorporating Novel Magnetism 15, 85 (2002). 
[28] Y. Zhao, W. T. Geng, K. T. Park and A. J. Freeman, 

Phys. Rev. B 64, 035207 (2001). 
[29] W. J. M. de Jonge and H. J. M. Swagten, J. Mag. Magn. 

Mat. 100, 322 (1991). 
[30] H. Asklund, L. liver, J. Kanski, J. Sadowski and R. Math- 

ieu, Phys. Rev. B 66, 115319 (2002). 
[31] J. Okabayashi, A. Kimura, T. Mizokawa, A. Fujimori, T. 

Hayashi and M. Tanaka, Phys. Rev. B 59, R2486 (1999). 
[32] J. H. Park, S. K. Kwon and B. L Min, Physica B281 & 

282, 703 (2000). 
[33] T. Dietl, F. Matsukara and H. Ohno, Phys. Rev. B 66, 

033203 (2002). 



